A SIMPLE SAMPLING METHOD FOR METRIC 
MEASURE SPACES 



EMIL SAUCAN 

Abstract. We introduce a new, simple metric method of sampling 
metric measure spaces, based on a well-known "snowflakeing operator" 
and we show that, as a consequence of a classical result of Assouad, the 
sampling of doubling metric spaces is bilipschitz equivalent to that of 
subsets of some R . Moreover, we compare this new method with two 
other approaches, in particular to one that represents a direct applica- 
tion of our triangulation method of metric measure spaces satisfying a 
generalized Ricci curvature condition. 



1. Introduction 

Sampling theory has its origins in the broad field of electrical engineering, 
more precisely in communication theory, namely in the works of Kotelnikov 
[30] . Nyquist [33] and mainly in the seminal papers of Shannon [65], [66] , 
[67] . Naturally, further developments along these lines followed, e.g. [32] . 
[38] , [69] , as well as natural applications to signal and image processing (see 
[72] for an extensive overview^), and also to the related field of graphics 
(e.g. [15] . [16] . [33]). Other applications include (but are not restricted to) 
information theory [73] and learning [69j. However, sampling theory rapidly 
transcended the boundaries of these applicative fields, to become a field of 
separate and sustained interest - see, e.g. [5T] . [68] . [35] . [36], amongst 
others. 
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Recently, driven mainly by what has become by now a common practice 
amongst the image processing community, namely to regard images as Rie- 
mannian manifolds embedded in higher dimensional spaces (usually M. n or 
§ n ) - see e.g. [22], [70], [63], [27], [55] - a geometric approach to sampling 
and its implications has emerged [26], [56], [57], [52] . 

A quite recent theoretical development concerns manifolds with densities 
[39] . [H], and even the more general metric measure spaces [2D] . [35] . [7T] . 
Such objects are not of pure mathematical interest, they arise naturally in a 
number of applicative fields. Amongst these applications, we first mention 
what is, perhaps, the the most immediate one, namely that of imaging. 
Indeed, classical grayscale (so called "natural") images can be viewed as 
densities over, say, the unit square. "Weighted" manifolds arise naturally 
in medical imaging, for the density of many types of MRI images is equal 
to the very proton density. Therefore, a correct modeling of such images, 
starting from the basic steps of sampling and reconstruction, will arguably 
produce more accurate results than the ones obtained with present methods. 
Manifold with densities are also encountered in information geometry [3], 
manifold learning, pattern recognition (both of these in conjunction with 
imaging) [33], [2], [32], bioinformatics [54J, graphics [28] and, perhaps most 
naturally, communication networks [53]. (In this last case, node and edge 
weights represent the relevant, characteristic features of the network.) 

The purpose of this article is rather straightforward: To introduce a sim- 
ple sampling method for metric measure spaces and to compare it to other 
approaches. Fittingly, the structure of the paper is also quite simple: In the 
following section we discuss two other sampling methods, the accent being 
placed upon the one stemming from our previous curvature-based triangu- 
lation method of metric measure spaces satisfying a generalized Ricci 
curvature condition, as introduced by Lott and Villani [5S] and Sturm |71| . 
In particular, we show that, indeed, the vertices of the said triangulation rep- 
resent a sampling of the given space. Moreover, based on the construction 
employed in [5T|, we also bring a result regarding the topological dimen- 
sion of weak CD(K,N) spaces. Section 3 represents the heart of paper: 
We introduce here the new sampling method, based upon of a well known 
"snowflaking operator" , and we show, as a consequence of a classical result 



SAMPLINGMMSP 



3 



of Assouad [I], [5], that sampling of doubling metric spaces is equivalent 
to that of subsets of R , for some N. We also prove that compact weak 
CD(K, N) spaces are Ahlfors A-regular, for K > 0. Section 4 represents 
the final section and is dedicated to a concluding discussion and comparison 
of the considered sampling methods. For the readers' convenience and for 
the sake of the paper's self containment, we bring some background material 
regarding curvature of metric measure spaces in an appendix. 

As it is clear from the very title and from the brief overview of the history 
of the subject given above, this paper is motivated not the least by the 
applicative goals, as they emerge in such fields as information geometry, 
image processing (and in particular, medical imaging), manifold learning, 
etc. Therefore, it is only natural that, while this paper is mathematical, not 
least in the origins of the suggested methods, we shall mention occasionally, 
in the relevant context, but mainly in the last section, a number of suggested 
applications. 

2. Other Approaches 

We present below two other methods of sampling metric measure spaces, 
besides the simple one announced in the title. We do not pretend that we 
thus exhaust all the possible approaches to the sampling problem on such 
manifolds, but rather concentrate on those methods that either are familiar 
to us or we deem to be more natural and/or important. 

2.1. Generalized Ricci Curvature-Based Sampling. Since the sam- 
pling method based upon the generalized curvature condition of Lott, Vil- 
lani and Sturm does not represent the one of main interest here, and since 
the technical definitions involved are rather lengthy and involved, we do 
not dwell upon them here, so not to disrupt the cursiveness of the exposi- 
tion. However, to preserve the self-containment of the paper, we adopt a 
compromise between fluency and clarity, by bringing them in an appendix. 

Considering that the exposition in [51] was rather roundabout and some- 
what didactic, and for the self containment of the paper, we bring here a 
concise proof. For this, we shall need the following notions: 

Definition 2.1. Let (X, d) be a metric space and let p%, . . . ,p no be points 
G X, satisfying the following conditions: 
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(1) The set {p\, . . . ,p no } is an e-net on X, i.e. the balls f3 n (pk,e), k = 
1, . . . , no cover X; 

(2) The balls /3 n (pk,e/2) are pairwise disjoint. 

Then the set {p±, . . . ,p no } is called a minimal e-net and the packing with 
the balls f3 n (pk,e/2), k = 1, . . . , no, is called an efficient packing. The set 
{(k, I) | k, I = 1, . . . ,no and f3 n (pi c ,e)r\(3 n (pi,e) ^ 0} is called the intersection 
pattern of the minimal e-net (of the efficient packing) . 

We our proof begin with the following lemmas: 

Lemma 2.2. Let (X, d, v) be a compact weak CD(K,N) space, N < oo, 
such that Supp^ = X and such that diamX < D. Then there exists n\ = 
n\(K, N, D), such that if {p±, . . . ,p no } is a minimal e-net in X , then uq < 
n\. 

Remark 2.3. Note that, since N < oo, the condition Suppz^ = X imposes 
no real restriction on X (see [73j, Theorem 30.2 and Remark 30.3). 

Lemma 2.4. Let (X,d,v) be a compact weak CD(K,N) space, N < oo, 
such that Suppz^ = X and such that diamX < D. Then there exists n<i = 
n2(N,K,D), such that, for any x G M n , \{j \j = 1, . . . ,no and /3 n (x,e) n 
(3 n (pj,e) 7^ 0}| < n%, for any minimal e-net {pi, . . . ,p no }- 

It is important to note that the integer n2 in the lemma above does not 
depend upon e. 

Lemma 2.5. Let {X\,d\,v\) and {X<i,di,v<i) be as in Lemma \2. 6 A and let 
{pi, . . . ,p no } and {qi, . . . , q no } be minimal e-nets with the same intersec- 
tion pattern, on X\, X<i, respectively. Then there exists a constant n% = 
ns(N, K, D,C), such that if d\{pi,pj) < C ■ e, then d2(qi,qj) < n^- e. 

The purpose of the lemmas above is to allow the construction of a trian- 
gulation of a weak (compact) metric measure space, as follows. Construct a 
simplicial complex having as vertices the centers of the balls /3 n (pk, e), in the 
following manner: Edges are connecting the centers of adjacent balls; fur- 
ther edges being added to ensure the cell complex obtained is triangulated 
to obtain a simplicial complex. (A concise and elegant exposition of these 
ideas in the classical (geometric differential) context, can be found in [9].) 
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The construction is possible since, by definition, weak CD(iT, N) spaces are 
geodesic. Moreover, in nonbranching spaces, the geodesies connecting two 
vertices of the triangulation are unique a.e. (see, e.g. [73], Theorem 30.17). 
(Recall that a geodesic metric space X is called nonbranching iff any two 
geodesies 71,72 : [0, i] — > X that coincide on a subinterval [0, rjo],0 < to < t, 
coincide on [0, t].) 

Remark 2.6. In smooth metric measure spaces (see the Appendix) one can 
actually produce a convex triangulation by choosing e to equal the con- 
vexity radius ConvRad(M") = inf{r > 0\f3 n (x,r) is convex, for all x £ 
M n }, and, moreover, control the convexity radius via the injectivity radius 
InjRad(M n ) = inf{Inj(x) | x G M n }, where Inj(x) = sup {r | &X-Vx\B n (x,r) ls a 
diffeomorphism}. What renders the construction above into a simple and 
practical triangulation method is a classical result of Cheeger, similar, later 
ones (see, e.g. |7J) showing that, in this case, there exists a universal posi- 
tive lower bound for InjRad(M) in terms of k, D and v, where v is the lower 
bound for the volume of M. 

However, many weak CD (if, N) spaces of interest fail to be locally convex, 
rendering the construction above impossible. Fortunately, local convexity 
does hold for an important class of metric measure spaces: Indeed, by [25], 
[75j . Alexfif] C CD((m — 1)K, m), where Alex [if] denotes the class of Tri- 
dimensional Alexandrov spaces with curvature > K (see [ID], [2DJ), equipped 
with the volume measure. 

We can formalize the construction above as: 

Theorem 2.7. (X, d, v) be a compact weak CD(K,N) space. Then X is 
triangulable. Moreover, if X is locally convex, in particular if it is a Tri- 
dimensional Alexandrov spaces with curvature > K , equipped with the vol- 
ume measure, one can ensure that the simplices of the triangulation are 
convex. 

Remark 2.8. For further improvements of the triangulation above and their 
applications, see [5lj . 

Of course, one is still has to ask himself if this triangulation result also 
represents, indeed, a sampling method? 
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It is quite classical by now that it represents a sampling method 
at least in the metric sense. More precisely, one assuredly obtains 
the Gromov-Hausdorff convergence of the e-net (viewed as a metric 
space) to M when e -> (see [20]. [ID]). 

Moreover, the same e-net also converges in measure, for instance in 
the measured Gromov-Hausdorff topology (see, e.g. [73]). For the 
choice of the weights in this context, see the next item. For a more 
detailed discussion regarding sampling based upon measures see the 
following subsection. 

However, given the fact that the vertices of the triangulation (points 
of the e-net) were constructed using generalized Ricci curvature, and 
since for the types of convergence above one hardly needs such so- 
phisticated tools, one is conducted to ask whether our method also 
assures convergence in the curvatures sense, i.e. that it is something 
akin - at least in spirit - to more classical results for PL approxima- 
tions of manifolds (such as those of [12] ) or of the Gromov-Hausdorff 
convergence for Alexandrov spaces (see [20], [ID]). 
The answer proves to be positive in this case, too, up to some (tech- 
nical) required adaptation of the notion of CT>(K, N) spaces to the 
"discrete" case of graphs (see [8j). First, let us define precisely the 
weights (i.e. the measure) for the considered space: Since the balls 
/3 n (pfc,e) cover M n , any sequence V(e m ), of efficient packings such 
that £ m — > when m — > oo, generates a discretization (X m ,d,v m ), 
in the sense of [8]: Consider the Dirichlet (Voronoi) cell complex 
(tesselation) <L m = {C(p m; k)k} of centers pk, m and atomic masses 
v m (p m ,k) = v[C(pm,k)]- Then taking X m = {p m ,k}, d the original 
metric of X and v m as defined, provides us with the said discretiza- 
tion. It follows, by |B], Theorem 4.1, that if Vol(M n ) < oo, the 
sequence (X m , d, v m ) converges in the W2 metric (see [71], [73]) to a 
metric measure space and, moreover, if (X, d, u) is a weak CD(K, N) 
space, then, for small enough e, so will be (X m , d, Um), but only in a 
generalized ("rough") sense. It should be noted that, by [8], Theorem 
3.10, the converse result also holds. More precisely, a (not necessar- 
ily compact) metric measure space (X, d, v) is a weak CD(K, N) 
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space if there exists a family {(X L ,d L ,i' L )}i of weak "discrete" met- 
ric measure CD(K L ,N L ) spaces, that converges to (X, d, v) (in the 
Wasserstein metric), where K b is in the weak (or "rough") sense, 
and diamX t < Dq, for some positive Dq, and such that K L — > K, 
when l —> 0. (For more details regarding the precise definition of 
rough curvature bounds and the proof of this and other related re- 
sults, see [8].) 

• The somewhat opposite problem is also of interest in some applica- 
tions, and in particular in image processing and graphics. Namely, if 
X is a topological manifold, such that each of the coordinate patches 
satisfies a weak CD(K L , N) condition, K L > Kq, one can triangulate 
each of these coordinate patches, and "glue" these local triangula- 
tions to obtain a global triangulation. The question is, naturally, 
whether the resulting triangulation will be a CD(Kq, N) space? The 
answer is positive, due to [6], Theorem 5.1, that shows that the local 
CD(K, N) condition implies a global, albeit slightly weaker one, and 
to [23], Theorem 4, where for a proof of the local-to-global property 
of the CD(K,N) condition is given. (A local-to-global property is 
also proved in [6] for a so called reduced CD(K,N) condition.) 

Before we bring the proofs of the lemmas, let us note that, incidentally, 
from Lemma 12.41 we obtain the following resulj§ (compare with |73j . Corol- 
lary 30.14.): 

Corollary 2.9. (X, d,u) be a compact weak CD(K,N) space. Then X has 
topological dimension < n2, where n2 = ri2(N,K,D) is as in Lemma \2.4\ 
above. 

Proof. The corollary follows immediately from the following characteriza- 
tion of topological dimension - see [25], Corollary to Theorem V 8. (p. 67) 
- namely that a compact space has dimension < n iff it has coverings of arbi- 
trarily small mesh and order < n. The definition (see, e.g. [25j, Definition V 
1.) of the order of a covering requires precisely the property guaranteed by 
Lemma [2.41 To show that compact weak CD(X, N) spaces enjoy the second 
required property, i.e. existence of coverings with arbitrarily small meshes, 
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one can proceed, amongst other possibilities, as follows: Start with the tri- 
angulation provided by Theorem 12,81 above. Since X is geodesic, one can 
built, iteratively, the barycentri<J§ subdivisions of any order, thus ensuring 
that the mesh of the triangulations tends to zero. Consider the balls having 
as diameters the sides of the obtained triangulation. Then, by a classical 
argument (see, for instance [7], [9]), this collection of balls represent, indeed, 
the required covering. □ 

The main (in fact, the only essential) tool in proving the Lemmas 2.1-2.3 
above is, by straightforward analogy with the classical case (see |21| ) the 
following generalized version of the Bishop-Gromov Comparison Theorem: 

Theorem 2.10 (Bishop-Gromov Inequality for Metric Measure Spaces, 
[7T]). Let (X,d,u) be a weak CD(K,N) space, N < oo, and let xq € SuppzA 
Then, for any r > 0, v{B{xQ,r)) = v(B[xQ,r\). Moreover, 

is a nonincreasing function of r, where 

N-l 



(2.2) S%{t) 



sin \/ fat) ifK>0 



t N ^ ifK = 



(sinh^Lt) ifK<0 



(Here, B(xo,r),B[xo,r] represent the standard convention for the open, re- 
spective closed ball center xq and radius r.) 

Proof of Lemma 2.1. Let {pi, . . . ,p no } be a minimal e-net on M n and let p 
be a point in MJ} - the A:-space form. Then, by Theorem 2.7 above 

v(B(p,r)) Volgfor) 

v(B(p,R)) ~ YolB(p,R) ' r ' 

for any p € M n . 

Let iq such that v(B(pi ,e/2)) is minimal. By 2.1.(2) it follows that 



^Since the "thickness" of the triangles (see, e.g. [51] for the definition) is not an issue 
here, the simple to produce barycentric subdivision will suffice. 
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v{M n ) YolB(p,D) 

n ° ~ u(B( Pl0 ,e/2)) ~ \o\B(p,e/2) ' 
(To obtain the last inequality, just take, in Bishop-Gromov Theorem, R 
diamikT < D.) 

The desired conclusion now follows by taking 



ni 



Vo\B(p,D) 
VolB(p, e/2) 



□ 



Proof of Lemma 2.3. Let j±, ...,j s be such that B(x,e)(lB(pj i ,e) ^0. Then 
B( Pji ,e/2) cB(x,5e/2). 

Let k 6 be such that B(pj k ,e/2) has minimal measure. Then 

(as in the proof of Lemma 2.2) it follows that: 

u(B(x,5e/2)) u(B( Pjk ,9e/2)) VolB(p,9e/2) 
S ~ u(B(p k ,e/2)) ~ v(B(p jk ,e/2)) ~ Vo\B(p,e/2) ' 
where p is as in the proof of the previous lemma. But 



VoL3(p,9e/2) _ ff' 2 S n K {r)d 



r 



VoW(p,s/2) $ £ J 2 S n K {r)dr ' 

and the function 

extends to a continuous function h : [0, D] — > R + , (since h{e) — > when 
e -> 0). 

□ 

Remark 2.11. It is important to note, both for the basic construction, and 
not least, for Corollary 12.91 that ri2 is independent of e. 

Proof of Lemma 2.4- Evidently, since d(pi,pj) < C ■ e, it follows that pj G 
B(pi,C ■ e). Thus, precisely as in the proof of the previous lemma, it follows 
that there exists n' = n'(C), 

, (r) Jo^ 2 Sl{r)dr 
n (C I = max — , 

fJ 2 Sl{r)dr 
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such that at most n' of the balls B(p±, e/2), B(p n r, e/2) are included in 
B(Pi, (C + i)e/2. 

Since {pi, . . . ,p no } and {qi, ■ ■ ■ , <?n } have the same intersection pattern, 
it follows that d(qi,qj) < n^{C), where ns(C) = 2[n'(C) — 1]. 

□ 

2.2. Gromov-Prohorov Metric-Based Sampling ("Measure Decides"). 

The approach above is natural, perhaps, from the point of view of a dif- 
ferential geometer since, as already mentioned, it extends to more general 
spaces and notions of curvature established ideas and techniques of classi- 
cal ("proper") differential geometry. It is however, an involved approach, 
and certainly not one that will appeal to people working in probability the- 
ory, statistic (or even main stream image processing, where histograms are 
a main "staple"). Such researchers will ask themselves, whether it is not 
possible to sample a metric measure space using solely the measure. The 
answer is, positive, as we shall see below, with the proviso that the metric 
should somehow be involved, since the goal is to find a sampling method 
for metric measure spaces. Therefore, the proper question would be if there 
exists such a method in which the main role is played by the measure (i.e. 
"measure decides"). 

We first have to define a "good" distance between metric measure spaces, 
where by "good" we mean here that it satisfies the requirement discussed 
above. Such a distance, would be, for instance, the Gromov-Prokhorov 
distance dap- Given two metric measure spaces X = (X,d,p) and y = 
(Y, p, u), we define 

(2.3) d GP (X,y) =infd P GuV)- 

Here the infimum is taken over all the measure preserving isometric embed- 
dings / : (X,d,p) — > (Z,5,X), g : (Y,p,u) — > (Z, 5, A), where (Z,S, A) is a 
common metric measure space (so the filiation of the Gromov-Prokhorov 
distance from the by now classical Gromov-Hausdorff metric - see below 
- is evident); and where dp denotes the classical Prokhorov distance, that 
can be defined in a very geometric manner (in the sense that is a straight 
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forward generalization of the Haussdorf metric - see (|2.6p below) as 
(2.4) 

d P (f/, u') = inf{r > | n'{F) < v'(N(F))+r , v'{F) < p'{M(F))+r , VF = F C 

where J\f(F) denotes the r-neighbourhood of F. 

This approach corresponds, according to [73J . to the "mainly measure" 
definition of isometry of metric measure spaces, namely that <p : (X, d, p) — > 
(Y, p, u) is an isometry of metric measure spaces iff it is a measure preserving 
isometry of the metric underlying spaces (Supp/z, d|s upP /J and (Suppi^, p\suppv)- 
Note that the metric indeed plays a role - albeit "subdued" - since the in- 
fimum is taken solely over isometric embeddings (and, moreover the metric 
does play a role in the "background" definition of the Prokhorov distance). 

Again, according to |73j . the "metric and measure, (but mainly metric)" 
approach is embodied in the Gromov-Hausdorff- Prokhorov distance dcHP, 
where 

(2.5) d GHP (X, y) = mf{d H (X, Y) + d P (X, y)} , 
where dn denotes the Hausdorff distance 

(2.6) d H (A, B) = inf{r > | A C B and B C A} , 

and all the other notations are as above, the infimum in (|2.5p being taken 
as in the definition of the Gromov-Prokhorov distance. 

In this setting, the fitting notion of isometry of metric measure spaces 
is that of isomorphism (or measure-preserving isometry) of metric measure 
spaces, meaning that, given (X, d, fi) and (Y, p, v), there exists a measurable 
bijection ip : X — > y such that (a) it is an isometry between (X, d) and 
(Y, p); and (b) it preserves measure, that is ip$p = v, where, as usual, VttA* 
denotes the push-forward of the measure p, i.e. ip$p(B) = p (■0 _1 (-B)), for 
any Borel set B C Y§ 

It turns out (see, for instance [73J, Theorem 27.26), that at least for 
the significant class of doubling spaces (see Definition 13.81 below), the two 
approaches are equivalent, and thus interchangeable in any mathematical or 
practical application. 



For other possible distances between metric measure spaces, see |73| . pp. 770-771. 
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3. "Snowflaking Operator" -Based Sampling ("Metric Decides") 

However simple and alluring the probabilistic approach may appear, to 
the geometer it seems somewhat unnatural. It is even less palatable to those 
whose interest is drove mainly by possible implementations, e.g. people 
working in information geometry, image processing, manifold learning, etc. 

Therefore, it is a natural desire to find a new metric that encapsulates 
the behaviors of both the original metric and of the given measure, at least 
as far as sampling (via e-nets) is concerned. 

3.1. Background: Quasimetrics and Doubling Measures. As general 
bibliographical references for the material in this subsection, including miss- 
ing proofs, we have used [23], [62], [63] . 

3.1.1. Quasimetrics. We begin with the following basic definition: 

Definition 3.1. Let X be a nonempty set. q : X x X — > M + is called a 
K-quasimetric iff 

(1) q(x,y) = iff x = 0; 

(2) q(x,y) = q(y,x), for any x,y G X; 

(3) q(x, y) < K(q(x, z) + q(z, y)), for any x, y, z G X . 

Remark 3.2. Some authors replace condition (2) above by the following 
weaker one: There exists Co > 1 such that q(x,y) < Coq(y,x), for any 
x,y G X. 

Remark 3.3. A number of brief comments: 

• A quasimetric is not necessarily a metric (while obviously, any metric 
is a quasimetric with K = 1). 

Counterexample 3.4. The following counterexample is not only the 
basic one, it is very important to us in the sequel: 

(3.1) q s (x,y) = (d(x,y)) s 

is a quasimetric for any s > 0, but not, in general, a metric, for 
s > 1. 

• Quasimetric balls can be defined precisely like metric balls, and the 
constitute the basis for a topology on X. 
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• For the next remark we need a definition that may appear a bit 
superfluous at this point, but it will prove to be highly relevant later 
on: 

Definition 3.5. Let (X, q) and (Y, p) be quasimetric spaces, and 
let / : X — > Y be an injection. / is called r]- quasisymmetric, where 
77 : [0, oo) — > [0, oo) is a homeomorphism iff 
, 3 2 ) p(/(x),/(a)) < ( q(x, a) 



p(f(x)J(b)) ~ '\q(x,b) 
for any distinct points x,a,b £ X. 

Intuitively, while quasisymmetric mappings may change the size of 
balls quite dramatically, they do not change very much their shape. 
This fact is important in the next proposition (see, e.g. [62J), that 
shows that whereas, as we noted above, q s is not a metric, the canon- 
ical injection (X,d) <-)■ (X,q s ) is quasisymmetric. 

Proposition 3.6. Let q be a K -quasimetric on X. Then, there 
exists so = sq(K) such that, for any < s < so there exists a metric 
d s on X, and a constant C = C(s,K) > 1, such that 

(3.3) —q s (x,y)<d s (x,y)<Cq s (x,y), 

where q s is as in \3.1\) , i.e. q s {x,y) = {q{x,y)) s . 

Remark 3.7. If q is a iT-quasimetric (K > 1), then q s is bilipschitz 
equivalent to d s , for any s > 0, such that (2K) 2s < 2, that is for any 
s > such that 

s< i(iog 2 tf + i). 

Moreover, the bilipschitz constant can be chosen to be 
C = (2K) 2s . 

3.1.2. From doubling measures to quasimetrics. We first remind the reader 
the following basic definition: 

Definition 3.8. Let (X, d, p) be a metric measure space X is called doubling 
iff [i is doubling, i.e. iff there exists a constant D such that, for any x £ X 
and any r > 0, 

(3.4) /i(fl d [x,2r]) <D(i(B d [x,r}) . 
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(Here Ba[x,r] denotes - as it standardly does - the closed ball of radius r, 
in the metric d.) A metric measure space (X,d,fx), where fi is doubling is 
sometimes called of homogenuous type. 

For the record, a metric measure space is a triple X = (X, d, fx) where 
(X, d) is a metric space and \x is a Borel measure on X. 

Remark 3.9. If (X,d,fx) is doubling, then it admits atoms only at isolated 
points. 

Remark 3.10. The connection between this definition and Section [2.11 be- 
yond the basic goal in both approaches, resides in the fact that any Rie- 
mannian manifold of nonnegative Ricci curvature is doubling (with respect 
to the volume measure) - see, e.g. [37] ■ (Indeed, it may be that this case 
represents one of the original motivations for studying doubling spaces.) 

Moreover, weak CD(K,N) spaces are locally doubling (on their support) 
- see [73J, Corollary 30.14 and globally doubling if X has bounded diameter. 
Since, by [73], Theorem 29.9, smooth weak CD(K, N) spaces are CD(K, N), 
the same assertions are true for smooth metric measure spaces. In fact, a 
stronger statement holds, since, by [73J, Corollary 18.11, smooth metric 
measure spaces, with 1 < N < oo, are globally doubling. 

Since this fact is not less important for our purposes, we remind the reader 
that if X is doubling (with added provisos of being Polish and compact) 
then, by [73J, Proposition 27.26 and Corollary 27.28, both approaches to 
convergence of metric measure spaces discussed in Section 2.2 are, indeed, 
equivalent, as far as convergence of e-nets - hence sampling - is concerned. 

For any s > 0, we define the quasimetric q^ tS as 



(This can be written in compact form as q^ jS (x,y) = (fx(B Xjy )) s , where 
B X)V = B[x, d(x, y)] U B[y, d(x, y)\.) 

Example 3.11. If X = M n , with fx = Vol n , and if s = 1/n, then q^ )S 
const • d Euc i- (In particular, for n = 2, q^ s = ^-d E ud ■) 



(3.5) 
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Remark 3.12. For X = W 1 , one can define q^ jS (x,y) simply by q )ltS (x,y) = 
(B[m, ^y^\)Y ' , where m denotes the midpoint of the segment xy. How- 
ever, in the general case, and in particular for graphs, one has to use the 
more general expression (|3.5p . 

Note that, if K is the quasimetric constant of q^s, then K = K(/i, s). 

Also, by Proposition 13.61 there exists sq = so(fi) > 0, such that q^ :S is 
bilipschitz equivalent to a metric cL S! for any < s < sq. This fact will 
play a crucial role in the remainder of the paper. 

Remark 3.13. Obviously, the geometry induced by the quasimetric q^g, and 
a fortiori by the metric cL s , will diverge widely from the geometry given by 
the original metric d. This is most evident in the properties of the "new" 
geodesies, in comparison with the "old" ones (e.g. when X = W 1 equipped 
with the standard Euclidean metric and with /i being the volume element.) 
However, the deformation of the geometry produced by (|3.5p is controlled, 
and many essential properties are preserved. (For further details, see [62], 
[63].) 

With the risk of being a bit confusing, but to be more specific, we hencefor- 
ward denote by D^^ the metric d^ s assured by Proposition ^. 61 for q s = q^ tS 
defined in (|3.5|) above. 

3.2. Equivalence of Nets in the Two Metrics (d and -D M . S ). The basic 
question for us is: 

Question 1. What is the relation between e-nets in d and in D^ s ? 

This can be decomposed into two more concrete questions: 
Question 2. Let Md be an d-e-net. Is it also a D^ tS -£-net? 
and the more interesting, for us 

Question 3. Is a D^g-e-net Md^ s also a d-e-net? More important, does it 
provide a sampling for [i as well? 

No general answer is available yet. However, we shall show that a positive 
answer exists for both questions in the important special case of Ahlfors 
regular spaces: 



16 



EMIL SAUCAN 



Definition 3.14. Let (X, d, f£) be a metric measure space, where [i is Borel 
regular. (X, d, /u) is called Ahlfors regular (of dimension a, or simply a- 
regular) iff there exists Co and a > such that 

(3.6) <fi(B[x,R])<C -R a , 
for any < R < diamX. 

Remark 3.15. Sometimes, the further hypothesis that (X,d) is complete is 
added, for convenience. 

Remark 3.16. For the a = 1/s, the geometry induced by the quasimetric 
Qfi,s, (or by the metric d^ s ) coincides, essentially, with the original geometry 
- see [BSJ , (see also Remark 13.131 above) . 

Imposing the quite mild Ahlfors regularity condition on a metric space 
assures that it "behaves in terms of size and mass distribution like Euclidean 
space"@. Moreover, compact Riemannian manifolds, endowed with their 
natural volume measure, are also Ahlfors regular (see, e.g., [UD], p. 273). 

Remark 3.17. It is immediate that an Ahlfors regular space is doubling. 

Remark 3.18. If (X,d,fi) is Ahlfors regular, then the set of atoms of [i 
is countable and, if x G A^, then there exists r > such that B(x, r) = {x}. 
(For a proof of this fact for, basically, Ahlfors 1-regular spaces, see [37] , 
Theorem 1.) 

Remark 3.19. It turns out that the specific measure fi in the definition above 
is not truly important, and in fact we can substitute for [i the Hausdorff 
measure H a . Indeed, for any Borel set E C X (and \i as above), there exists 
C > 1 such that 

(3.7) ^j-H a < fi(E) <C'-H a . 
Moreover, the Hausdorff dimension dim#(X) = a% 

5 [63], p. 15 

^For further results relating fj,, -D^.c and an arbitrary Ahlfors regular measure v on X, 
see [63], p. 57 ff. 
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The following result is, perhaps, of no great importance by itself, it is, 
however, significant in our context, since it relates between the first and 
third sampling techniques by demarking a class of spaces for which both of 
the mentioned techniques are applicable: 

Lemma 3.20. Let X = (X,d,(i) be a weak CD(K, N) compact space, K > 
0. Then X is Ahlfors N -regular. 

Proof. By [73], Corollary 30.12, that there exists C = C(K,N,R) such that 
the following holds 

n(B(x,r))>C-n(B(y,R))r N , 

for any x,y,r,R such that B{x,r) C B(y,R). In particular, for x = y and 
R = diampT), we obtain that fi(B(x,r)) > C ■ ^(X)r N , that is 

KB(x,r)) >C x r N , 

where C x = C ■ \i{X). 
The opposite inequality 

li(B(x,r)) < C 2 r N 

in (13. 6p follows immediately from the Bishop-Gromov inequality (12. 1QH for 
K > 0, where the constant C 2 is the constant that appears in the formula 
for the volume of the ball of radius r in the space form of dimension N, 
given as a a function of J* Q r Sj[(t)dt. 

One can obtain the precise form of the double inequality in (|3.6p by choos- 
ing a constant C* as C* = max{C2, Cj~ 1 }. □ 

We next prove now that in Ahlfors regular spaces, d-nets and D M)S -nets 
are, indeed, equivalent. More precisely, we can formulate the following 

Proposition 3.21. Let (X,d,fi) be a 1/s-regular space. Then d-nets and 
D^^-nets are equivalent. 
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Proof. In fact, we can prove the proposition above not only for the metric 
.D^s, but also for the more lax quasimetric q^^. Then assertion for the metric 
D^s follows from Proposition 13.61 (with different constants, of course). 



d-e-nets =>• g M)S -e:-nets This implication is quite easy and it holds for 
general a (not necessarily equal to l/s)0 Indeed, if d[x,y) < e\, for 
some £1 > 0, then q^ s (x,y) = (p,(B[x,d(x,y)}) + \i (B[y, d{x, y)]) ) s < 
{n{B[x,ex}) + /i( J B[y,ei]))" < C*e{ a , where C* = 2 a C%. 



q^-e-nets <= d-e-nets We begin by trying to better understand (following 
62] ) the geometry of g^-bails, which we denote by /3: /3[x, p) = {y | ?^ j£ (a?, y) < 
r}. (If p[x,p] = X, let p = inf{p' | /3[x,p] = X}.) Let r = inf{r' | 0[x,p] C 
£>[x, r']}. Then there exists z G /3[x, p] such that d(x, z) > r/2. It can be 
shown that there exists C3, independent of z and such that 

(3.8) B{x,r/C 3 )cp[x,p])CB[x,r]. 
Moreover, there exists C4 such that 

(3.9) ^r<p(B[x,r])- s <C 4 p. 

Using (|3.8p . (j3.9|) and the doubling condition (for we infer that 

(3.10) -?r<mx,p})- s <C iP . 

Now the remainder of the proof is quite elementary. Indeed, lUj) holds 
for y instead of x, therefore, for s > 1 we have: 

^r) S <KP[x,p])<(C4P)-\(£-) 8 <p(P[y,p])<(C 4P )- s ; 



hence 



thence 



that is 



_P_ 



2 ( £ ) < M(/3[a:,p]) + < 2(C 4 p)" 



< 



( / x(/3[x, /3 ])+ / i(/3[y,p])) S <2 s (C 4 p)- 



2 s (^Q <g M ,s<2 s (C4/>)- 



'''in fact, it holds for so called normal (metric measure) spaces, i.e. such that there exist 
< ci,C2 < 00, satisfying c\r < fi(B[x,r]) < C2T2, for all x £ X and for any r > 0, such 
that n({x}) < r < n{X). 
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For < s < 1, the same argument applies, because, while the inequalities 
change orientation, the double inequality of type (13. 6p still holds. 

□ 

3.3. Qui Prodest? The unavoidable question which we are confronted with 
is whether the sampling result above has any practical potential. This ques- 
tion rises not least because the notion of Ahlfors regular spaces it is perhaps 
less known to the sampling community, therefore it is natural to ask whether 
such spaces are fairly common or just yet another technical artifice. 

We shall show that, in fact, Ahfors regular spaces are quite abundant, 
and also how one can construct such spaces - in a manner that will appear, 
by now, quite natural. 

We first have, however, to introduce yet another definition, that ensures 
that a metric space (X, d) contains no "isolated islands" : 

Definition 3.22. A metric space (X,d) is called uniformly perfect iff there 
exists C\ > such that, for any x £ X and any < r < diamX, there exists 
y € X such that 

(3.11) -^<d(x,y)<r. 
(For other, equivalent definitions, see [23J, |62j.) 

One can show that, while the basic doubling condition ensures that balls 
do not grow too fast, in uniformly perfect spaces they also do not decrease 
at a too steep rate, more precisely that 

(3.12) p(B[x, a k r}) < (1 - a) h fi(B[x, r]) , 

for any k € N, any x € X and any < r < diamX. 

Using this fact, it is easy to prove (see [62], [23]) the following result, that 
shows how to canonically construct Ahlfors regular spaces: 

Proposition 3.23. Let (X, d) a uniformly perfect metric space and let [i 
be a doubling measure on X. Then (X, d^ s , //) is --Ahlfors regular for any 
s > 0. Moreover, there exists so > 0, such that the canonical injection 
(X, d) (X, cL s ) i> s quasisymmetric for any < s < sp- 
in particular, graphs endowed with metrics are --Ahlfors regular (for 
any s > 0), hence e-ci-nets in such graphs are also e-d^^-nets. 
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Remark 3.24. The proof of the first assertion of the proposition above is 
basically included in the second part of the proof of Proposition 13.211 

Remark 3.25. As a consequence of (13. 12ft . it follows that doubling measures 
on uniformly perfect spaces have no atoms. (Compare to Remarks 13.91 and 
CHS above.) 

In the opposite direction, it is also easy to prove that Ahlfors regular 
spaces are uniformly perfect (see [62]). In fact, we have the following result 
(for a proof, see [62]): 

Proposition 3.26. A metric space (X, d) is quasisymmetrically equivalent 
to an Ahlfors regular space iff (X, d) is uniformly perfect and admits a dou- 
bling measure. 

So Ahfors regular spaces represent controlled (quasisymmetric) deforma- 
tions of "nice" doubling spaces, where "nicety" is formulated in terms of 
the rather mild restriction given by uniform perfectness. However, the sec- 
ond condition in Proposition 13.261 above, on the existence of the doubling 
measure, may be difficult to verify. Fortunately, a simpler, purely metric 
condition is equivalent to it, at least for complete spaces, namely: 

Definition 3.27. A metric space (X, d) is called doubling iff there exists 
D\ > 1, such that any ball in X, of radius r, can be covered by at most D\ 
balls of radius r/2. 

(Obviously, there is nothing special about balls, and the metric doubling 
condition can be formulated in terms of general sets of bounded diameter.) 

Not surprisingly, there exists a connection between the notions of doubling 
metric and doubling measure. More precisely, we have the following 

Lemma 3.28. Let (X,d) be a metric space such that there exists a doubling 
measure fj, on X. Then (X,d) is doubling (as a metric space). 

(For a proof, see [23], [62].) 

The converse statement does not always true, a counterexample being 
provided by (Q, d euc i) (see [23], p. 103). However, it does hold for complete 
spaces: 
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Theorem 3.29 (Luukkainen-Saksman [36]). Let (X,d) be a doubling, com- 
plete metric space. Then X carries a doubling measure. 

Corollary 3.30. Any compact, doubling metric spaces carries a doubling 
measure \i. 

The corollary above obviously holds for finite graphs. 

Remark 3.31. If (X, d) is a-Ahlfors regular, then, by Remark 13.191 above, fi 
can be taken as the Hausdorff a-measure. 

From the above noted equivalence between complete doubling metrics and 
doubling measures, and from Proposition 13.261 we obtain the following 

Corollary 3.32. Any complete, uniformly perfect metric spaces is qua- 
sisymmetrically equivalent to a Ahlfors regular space. 

Before we formulate the important theorem of Assouad on which we shall 
base our sampling result, we give, for convenience, the following definition: 

Definition 3.33. If (X, d) is a metric space, then the metric space (X, d £ ),0 < 
e < 1, is called a snowflaked version of (X, d). 

Theorem 3.34 (Assouad, [I], [5]). Let (X,d) be a doubling metric space. 
Then, for each < e < 1, there exists N, such that its e-snowflaked version 
is bilipschitz equivalent to a subset ofM, N , quantitatively. 

Here, quantitatively means that the embedding dimension N and the bilip- 
schitz constant L depend solely on the doubling constant D of X and on 
the "snowflaking" factor e, i.e. 

N = N{D,e), L = L(D,e). 

Remark 3.35. Assouad's result does not hold, in general, for e = 1. (For a 
counterexample, see [23], p. 99). 

Combining Assouad's Theorem, Corollary 13.321 and our own sampling 
result, we can now enunciate the following sampling "meta-theorem" : 

Theorem 3.36. Sampling of Ahlfors regular metric measure spaces is qua- 
sisymmetrically equivalent, quantitatively, to the sampling of sets in~R N , for 
some N. 



22 



EMIL SAUCAN 



Remark 3.37. The beauty of Assouad's Theorem - and even more so its 
applicability in the sampling of real data - is marred by the "course of 
dimensionality": Given that N = N(D,e), the fear exist that, as in the 
case of Nash's Embedding Theorem [II], [32], the embedding dimension is 
prohibitively high for general manifolds (i.e. data). Obviously, this is even 
more important if low distortion - i.e. (bi-)lipschitz constant - is an im- 
perative (as it usually is), that is for e close to 0. And, indeed, Assouad's 
original construction provides lim e ^o N(D, e) = oo. So it would seem that, 
the price to pay for low distortion is a high embedding dimension. It is a 
quite recent result of Naor and Naiman [30] (itself based on ideas of Abra- 
ham, Bartal and Neiman [I]), that, in fact, given a (separable) D-doubling 
metric space, there exist N = N(D) € N and L = L(D,e), such that for 
any e £ (0,1/2), the (1 — e)-snowflaked version of X admits a bilipschitz 
embedding in M. N , with distortion L. Moreover, specific upper bounds for N 
and L are given: N < alogD,L < b ^ j , where a and b are constants. 
So it appears that, at least as far as Assouad's Theorem is concerned, the 
snowflaking-based embedding is feasible. 

At this point, one has to ask oneself whether this result can be improved. 
The belief in the possibility of such an improvement rests upon the following 
two facts: One one hand, Assouad's Theorem assures the existence of a 
bilipschitz embedding, which represents a much stronger condition then mere 
qusysymmetrjj^. On the other hand, as we have seen, Ahlfors rigidity is not 
the most easy property to check directly on a metric measure space, therefore 
one naturally would wish to find a sampling result similar to Theorem 13.361 
that would hold for general doubling spaces. Such a result does exist, and 
it makes appeal again to the quasimetric q^g as defined by (3.3). However, 
we have to make an additional assumption, that ensures that g^-lengths of 
curves in M. N do not "shrink" too much, due to the presence of the measure 
H in the definition of g^ iS (see [62]). We encode this restriction via 



However, quasisymmetry represents a much more flexible analytic tool, than the rigid 
bilipschitz condition - see [23], [62], [63] for a deeper and far more detailed discussion. 
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Definition 3.38. A doubling measure \x on M N is called a metric doubling 
measure iff there exist a constant Cq, and a metric 5, such that 

7^0, V) < Q u I < Ce^, y) , 
for any x,y £ X, where g i is associated to fi as in (3.3), with s = 1/n. 

We can now formulate the desired result, in terms of metric doubling 
measures: 

Theorem 3.39 (Semmes [59] Theorem 1.15, [62], Proposition B. 20.2). Let 

(X, d) be a doubling metric space. Then there exists a natural number N 
and a metric doubling measure u, such that (X, d) is bilipschitz equivalent 
to a subset of(M. ,q,, i ), where q, i is as above. 

This is a most encouraging result, and the idea of the proof is quite simple: 
By Assouad's Theorem, (M, d^) is bilipschitz equivalent to a subset Y of 
some W N . The sought for measure on W N will be define as \x = dist(x, Y n )dx 
- for details of the proof see [59J . 

One would naturally would hope that (R n , q i) can be bilipschitzly em- 
bedded in some R , for any doubling measure u. This is a quite ambitious 
wish and, unfortunately, it is not true in general (see |59j). However, such 
an embedding exists for "most" metric doubling measures - for a precis 
formulation and the proof see [59]. Still, we can formulate the fitting sam- 
pling result (recall that given the quasimetric g„ )S , there exists a metric d s 
bilipschitz equivalent to it): 

Theorem 3.40. Sampling of doubling metric spaces is bilipschitz equiva- 
lent quantitatively to the sampling of sets in (M. N ,d±), for some N, where 

N 

dj_ represents the snowflaked version of d, associated to a certain metric 

N 

doubling measure u. 

4. Discussion and Final Comments 

A number of concluding remarks, regarding the relative advantages of the 
three sampling methods exposed above, are mandatory. 

As far as simplicity is concerned, then obviously the third method is 
the preferred one, as emphasized already (even in the title of the article): 
It is the most intuitive (at least for a geometer), employing just a quite 
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simple metric. Of course, there exists a trade-off between precision and 
simplicity, ensuing from the fact that, if one insists on working with an 
actual metric and not a "mere" quasimetric, then he/she has to be content 
with approximation provided by Proposition 13.61 Moreover, it is highly 
adaptable, via the parameter s, that allows for sampling at different scale 
(as envisioned originally by Semmes, for quite different ends). This is not a 
negligible advantage, since in many application it is not a priori clear at what 
scale the data should be sampled at. Imaging data is the first example that 
cames to mind, see, e.g. [T7], [IB], [2S], [50]. Astronomy (in its cosmological 
setting) represents, probably another such case. This problem appears even 
more poignant for data where little information regarding the structure of 
the data exists, and in particular no natural dimension is available. Data 
from bioinformatics appertains to this category. While the second ( "measure 
decides") method is also relatively simple, it shares with the last approach 
one common weakness, that is the fact that, in the absence of of curvature, 
there exists no way of determining the metric density of the sampling points, 
akin to that of [11] (see also [56] for an application in imaging). Therefore, at 
this stage, this methods seem to be feasible only for data that is intrinsically 
"almost flat" (such as it appears in [15]). 

This highlights the relative advantage of the curvature-based method 
upon the other two: While it is, admittedly, conceptually the more com- 
plicate, and computationally quite involved, it is the only one that satisfies 
the density condition above and, more important, it is the only algorithmic 
one. 

However, on a more theoretical level, the superiority belongs to the metric 
method, since, due to Theorem 13.361 it reduces - even quantitatively - the 
sampling of quite general metric measure spaces to that of subsets of M. N . 
This is also important in applications, since, by employing such methods 
as those devised e.g. in [56J, more general kind of data, and not just the 
"almost flat" one can be sampled (albeit by making appeal to the extrinsic 
curvature of a specific embedding). 
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Appendix - Generalized Ricci Curvature of Metric Measure 

Spaces 

We bring here only the minimal amount of definitions needed as a back- 
ground material for Section 2.1. (The interested reader can consult, for 
further details, the exhaustive monograph |73j and, of course, the original 
papers [35] and [7T].) 



4.1. Smooth Metric Measure Spaces. Let M = M n be a complete, con- 
nected n-dimensional Riemannian manifold. One wishes to extend results 
regarding Ricci curvature to the case when M n is equipped with a measure 
that is not dVol. Usually (at least in our context) such a measure is taken 
to be of the form 

(4.1) u(dx) = e- y{x) Yo\{dx) , 

where V : M n ->• R, V G C 2 (R). Note also that any smooth positive 
probability measure can be written in this manner. Then (M,d,u), where 
d is the geodesic distance, is a metric measure space. 

Remark 4.1. A standard measure u, in the context of image processing (but 
not only) is the gaussian measure on M n : 

(4 2) >) = e ~ lX]2dx 

To preserve geometric significance of the Ricci tensor, one has to modify 
its definition as follows: 

(4.3 Ric Nu = Ric + V 2 V 

N — n 

Here VV<8> W is a quadratic form on TM n , and V 2 V is the Hessian matrix 
Hess, defined as: 

(4.4) (X7V ®X7V) x (v) = (VV{x) -v) 2 . 
Therefore 

(4.5) Ri Cjv (7) = (Ric + V 2 V)( 7 ) - (V J' 7)2 • 

TV — n 

Here N is the so called effective dimension and is to be inputed. 
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Remark 4.2. (i) If N < n then Ricjv,i/ = — oo 

(ii) If N = n then, by convention, x oo = 0, therefore (I4.3P is 
still defined even if W = 0, in particular Ric nj voi = Ric (since, in this case 
V = 0). 

(hi) If N = oo then Ric^ = Ric + V 2 V. 

The Ricci curvature boundedness condition of the classical Bishop-Gromov 
is paralleled in the case of smooth metric measure spaces by the following 
immediate generalization of the classical definition: 

Definition 4.3. (M, d, v) satisfies the curvature- dimension estimate CD(K, N) 
iff there exist K G M and N € [1, oo], such that Ric jsj v > K and n < N. (If 
v = dVol, then the first condition reduces to the classical Ric > K.) 

Remark AA. Intuitively, "M has dimension n but pretends to have dimen- 
sion N. (Identity theft)"§ 

The need for such a parametric dimension stems, in particular, from the 
desire to extend the Bishop-Gromov Theorem to metric spaces (or more 
precisely, to length spaces), for which no innate notion of dimension exists. 

Remark 4.5. For a number of equivalent conditions, see [73], Theorem 14.8. 

4.2. Weak CD(K,N) Spaces. 

Definition 4.6. Let (X, fi) and (Y, v) be two measure spaces. A coupling 
(or transference (transport) plan) of [i and v is a measure n on X x Y 
with marginals [i and v (on X and Y, respectively), i.e. such that, for all 
measurable sets A C X and B C Y, the following hold: n[A X Y] = fi[A] 
and ir[X X B] = v[B]. 

Definition 4.7. Let (X,/j>) and (Y,v) be as above and let c = c(x,y) be a 
(positive) cost function on X x Y. Consider the Monge-Kantorovich mini- 
mization problem: 

(4.6) inf / c(x,y)dir(x,y) , 

where the infimum is taken over all the transport plans. The transport plans 
attaining the infimum are called optimal transport (transference) plans. 

9 J. Lott IM]. 
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Before we can proceed, we must recall the following definition and facts: 

Definition 4.8. Let (X, d) be a Polish space, and let P(X) denote the set of 
Borel probability measures on X. Then the Wasserstein distance (of order 
2) on P(X) is defined as 



where the infimum is taken over all the transference plans between \x and v. 
Definition 4.9. The Wasserstein space is defined as 



where xo £ X is an arbitrary point. 

Remark 4.10. The definition above does not depend upon the choice of xq 
and W2 is a metric on P^{X). Moreover, if X is Polish (compact), P^{X) is 
also Polish (compact). 

Definition 4.11. Let (X,d) be a compact, geodesic, Polish space and let 
T = {7 : [0,1] — > X \ 7 a minimal geodesic}, and denote by et ■ T — >• X 
the (continuous) evaluation map, et : (7) = j(t). Let E : T — > X x X be 
defined as £"(7) = (60(7)161(7)). A dynamical transference plan is a pair 
(it, II), where ir is a transference plan and LI is a Borel measure, such that 
E#U = ir. (ir, LT) is called optimal if ir is optimal. 

Definition 4.12. Let IT be an optimal dynamical transference plan. Then 
the one-parameter family {t-H}te[o,i\> fH = ( e *)#n is called a displacement 
interpolation 

We can now quote the following result ([73], Theorem 7.21 and Corollary 
7.22), connecting the geometry of the Wasserstein space to classical mass 
transport: 

Proposition 4.13. Any displacement interpolation is a Wasserstein geo- 
desic, and conversely, any Wasserstein geodesic is obtained as a displace- 
ment interpolation from an optimal displacement interpolation. 



(4.7) 




(4.8) 
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Definition 4.14. Given N € [l,oo], the displacement convexity class T>Cn 
is defined as the set of convex, continuous functions U : M+ — >■ R, U € 
C 2 (M + \ {0}), such that [7(0) = and such that 

V*'*> r X-X/N 
is nondecreasing (as a function of r). 

Remark 4.15. For equivalent defining conditions for the class T>Cn see |73j . 
Definition 17.1. 



Definition 4.16. Let (X, d, v) be a a locally compact metric measure space, 
such that the measure v is locally finite, and let U be a continuous, convex 
function t7 : M+ -> M, C/ € C 2 (M+ \ {0}), such that U(0) = 0. Consider 
a measure /i on X, having compact support, and let \i = pv + \i s be its 
Lebesgue decomposition into absolutely continuous and singular parts. 

Then we define the (integral) functional U„ (with nonlinearity U and 
reference measure v) by 

(4.10) U v = f U(p(x))u(dx) + U'(oo)[i s [X] . 

Jx 

Moreover, if {7r(dy\x)} x ^x is a family of probability measures on X and 
if (3 : U x U — > (0, oo] is a measurable function, we define an (integral) 
functional Un,u (with nonlinearity U, reference measure v, coupling ir and 
distortion coefficient j3) by: 

(4.11) [/£ = / U (-^-) P(x,y)7r(dy\x)u(dx) + U'(oo)p JS [X}. 

Juxu \P(x,y)J 

Usually (e.g. in the definition of weak CD(K, N) spaces) (3 is taken to be 
the reference distortion coefficients: 

Definition 4.17. Let x,y be two points in a metric space (X,d), and con- 
sider the numbers K €, N £ [l,oo] and t £ [0, 1]. We define the reference 
distortion coefficients /3^ K ' N \x,y) as follows: 
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(1) If t E (0, 1] and 1< N < oo, then 

+oo if K > and a > 7r , 



(4.12) /3 f' JV )(x,y) = <' 



sin (tct) 
t sin a 



7V-1 



if K > and a E [0, it] 



if K = . 



where 



(4.13) 



a 



\K\ 



-d(x,y) 



N-1 

2) In the limit cases N — > 1 and iV — > oo, define 



(4.14) 

and 

(4.15) 

(3) If i = 0, then 
(4.16) 



A K ' l) (x,y) 



+oo if K > . 
1 if K < : 



/3f^)( X ,y) = e f(l"* 2 M^)- 



/3f' Af) (x, y ) = l 



Remark 4.18. If X is the model space for CB(K,N) (see [73] p. 387), then 
fi( K > N ) ig the distortion coefficient on X. 

We can now bring the definition we are interested in: 

Definition 4.19. Let (X, d, v) be a locally compact, complete, cr-finite 
metric measure geodesic space, and let K E M, iV E [l,oo]. We say that 
(X, d, i/) satisfies a weak CD(K, N) condition (or that it is a weak CD (if, N) 
space) iff for any two probability measures /io>Ml with compact supports 
Supp^i, Supp/i2 C Suppz/, there exist a displacement interpolation fJ<to<t<l 
and an associated optimal coupling ir of [Ao,m such that, for all C7 E PC at, 
and for all t E [0, 1], the following holds: 



(4.17) 



R (K,N) 
TP\—t 



U v (jH)<(l-t)UZy {/*>) + tU% v (mi 



(K,N) 
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(Here we denote tt = S#n, where S(x,y) = (y,x).) 

Remark 4.20. In fact, the geodesicity condition is somewhat superfluous, 
since a locally compact, complete metric space is geodesic (see, e.g. [23] , 
9.14). 
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